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Using the mimetic matter approach, we study F{R) gravity with scalar potential and Lagrange 
multiplier constraint. As we demonstrate, for a given F{R) gravity and for suitably chosen mimetic 
potential, it is possible to realize inflationary cosmology consistent with Planck observations. We 
also investigate the de Sitter solutions of the mimetic F{R) theory and study the stability of the 
solutions, when these exist, towards linear perturbations, with the unstable solutions, which can 
provide a mechanism for graceful exit from inflation. Finally, we describe a reconstruction method 
which can yield the F{R) gravity that can generate realistic inflationary cosmological evolution, 
given the mimetic potential and the Hubble rate. 
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I. INTRODUCTION 

Our Universe’s current evolution status is determined by the contribution of ordinary matter with 0^ 4.9%, 

cold dark matter with UIdm ^ 26.8% and also dark energy with UIde ^ 68.3%, with the dark energy contribution 
controlling the late-time acceleration, verified observationally in the late 90’s [l|. A consistent description of late-time 
acceleration is provided by the modified theories of gravity (see reviews, i)- There is number of modified theories 
of gravity that can model quite successfully dark energy (for discussion on dark energy, seeQ). On the contrary, 
for cold dark matter very few are known, apart from the fact that this form of matter modifies the galactic curves. 
In the literature, it is common to consider dark matter as a particle that has no-interaction with luminous matter. 
In the context of supersymmetric modifications of the Standard Model, one particle that is a potential candidate 
for describing cold dark matter, is the neutralino Q, but there exist also other particles that could describe dark 
matter, see for example @ and references therein. Apart from the particle nature of cold dark matter, there exists 
a recent alternative description to cold dark matte:^ provided by the so-called mimetic modification of General 
Relativity, proposed by Chamseddine and Mukhanov @ and further studied in 0,0 . In their paper 0] , Chamseddine 
and Mukhanov isolated the conformal degree of freedom of Einstein-Hilbert gravity in a covariant way, and in the 
resulting theory, the physical metric is defined with the account of an auxiliary scalar field, which appears through 
its first derivatives. In effect, the conformal degree of freedom becomes dynamical, even in the absence of any matter 
fluids, and more importantly this conformal degree of freedom can mimic cold dark matter. A natural extension of 
the mimetic Einstein-Hilbert gravity is provided by mimetic F{R) gravity 0], since these theories could potentially 
have the attributes of the mimetic dark matter approach, but it is possible to also describe early-time and late-time 
acceleration within the same theoretical framework [ 13 - 

In this letter we shall study Jordan frame mimetic F(R) gravity, with the mimetic dark matter part having also 
a potential term. For some informative reviews on F{R) theories we refer to 0 ], and for some early but worthy to 
mention studies, see m- We shall use the Lagrange multiplier approach [10 to modified gravity, so that we achieve a 
consistent incorporation of the mimetic matter potential, in the Jordan frame F{R) framework. Our aim is to provide 
a reconstruction technique for obtaining the mimetic potential which generates a cosmological evolution compatible 
with the recent Planck data [13, by using viable F(i?) theory, with viable meaning that they satisfy the local and 
large scale constraints posed on these by observations 0. We shall mainly work in the Jordan frame and in order to 
have a better command on the cosmological evolution we shall use the e-folding number N as a variable, and we shall 
express all the involved to the calculation quantities in terms of this variable. Also, we shall briefly provide a second 
reconstruction technique, which enables us to find which F(R) gravity can generate a specific cosmological evolution 
with a specific mimetic potential. Finally, we investigate the existence of unstable de Sitter solutions for some F{R) 
gravities with mimetic potential. 
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II. CONSISTENT COSMOLOGIES FROM LAGRANGE-MULTIPLIER MIMETIC F{R) GRAVITY 

In the context of mimetic approach, the conformal symmetry is considered as an internal degree of freedom and it 
is not violated. The physical metric 5 ^ 1 ^, is expressed in terms of an auxiliary metric and of an auxiliary scalar 
field (j), in the following way: 


9i^i' = ■ ( 1 ) 

Instead of considering the variation of the gravitational action with respect to the physical metric one can consider 
variation with respect to the auxiliary metric and with respect to </>. From Eq. ([T]) it easily follows that, 

(I>)dp4>d^(t> = - 1 . ( 2 ) 

The parametrization © is invariant under the Weyl transformation g^i, = and the auxiliary metric g^v does 

not appear explicitly in the action. We shall assume that the physical metric g^n is a flat Friedmann-Robertson-Walker 
(FRW) metric with line element. 


ds^ =—dt^ + a{tY ^ (da;*)^ , (3) 

i=l,2,3 

where the parameter a{t) denotes the scale factor, and the scalar (f) depends only on time. For the flat FRW background, 
the Ricci scalar has the form R = 6 (^H + 2H^^, where H{t) is the Hubble rate. The Jordan frame gravitational 
action for the mimetic F{R) theory with scalar potential V{(j)) and Lagrange multiplier A(^), is the following, 

S = J dx^y/^ (F {R{gp,)) - V{^) + A {g^’^dpcj^d,^ + 1 )) , (4) 

where for simplicity, no extra matter is present. By varying the action o, with respect to the physical metric g^u, 
we can easily obtain, 

^gp.F{R) - Rp,F'{R) + - gp,UF'{R) (5) 

^ {9^''dp(j)da(l) + 1)) - \dp(j)du4i = 0 . 

Varying the action of Eq. with respect to the auxiliary scalar (j), we obtain, 

-2V'^{Xdpcl>)-V'{cj)) = 0, ( 6 ) 

where the prime denotes differentiation with respect to the auxiliary scalar (j). Upon varying the action (jd]) with 
respect to A, we get. 


gP'^dp(t)d^(j)= -I, (7) 

which is identical to Eq. @. Assuming the flat FRW background of Eq. ©, and also that the scalar field depends 
only on the cosmic time t, the equations of motion ([5]), dH) and ([71) are written as follows, 

- F{R) + 6{H + H^)F'{R) - - A ((^2 ^ ^ ^ g , ( 8 ) 

at 

F{R) - 2{H + iH^) + - 1) - = 0 ’ (9) 

at^ at 


+ 6HX<P - V\(l)) = 0 , ( 10 ) 

at 


<^ 2-1 = 0 , 


( 11 ) 
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where the “dot” in all the above equations denotes differentiation with respect to the cosmic time t, while the “prime” 
in Eqs. ® and © denotes differentiation with respect to the scalar curvature, while in Eq. (HHD denotes differentiation 
with respect to the scalar field (j). Erom Eq. m, it follows that the scalar field (f) can be identified with the cosmic 
time t (the same identification was possible in the Einstein-Hilbert mimetic gravity, as an implication of Eq. ([5]) , see 
0). Since t = (j), Eq. can be written in the following way, 

F{R) - 2{H + 3H^) + - V{t) = 0 , (12) 

aF at 

from which it easily follows that the scalar field potential V{(j) = t) that generates the Hubble rate H(t) , for a specific 
F{R) gravity, is equal to. 


E((^ = t) = 2^-^^ + + F{R) - 2{H + 3H^). (13) 

dt^ dt 

Upon combining Eqs. m 3,nci (ISp . w6 obtnin the ctncilytic form of tliG LcigrcingG multipliGr, cis 3^ function of the cosmic 
time which reads, 


X{t) = -3H^^^^ + 3{H + H‘^)-^F{R). (14) 

Then, we have at hand a reconstruction method with which we can find which scalar potential and Lagrange multiplier 
can generate a specific Hubble rate, for a given arbitrarily chosen F{R) gravity. 

The novelty of this reconstruction method we propose is that, in principle we can choose a viable F{R) gravity 
and a Hubble rate that is compatible with the Planck data and find the corresponding mimetic scalar potential that 
generates the evolution, for the given F{R) gravity. In order to exemplify our method, we shall use an illustrative 
example by choosing an arbitrary F{R) gravity. Consider for example the Jordan frame cubic F{R) gravity, 

F{R) = R + dR^ , (15) 


with arbitrary parameter d. Our aim is to find which scalar field potential V{(j} = t) reproduces a cosmology compatible 
with the observations, for the F(R) gravity chosen as in Eq. (fTK)) . We shall use the slow-roll parameters and 
observational indices found in Ref. [l4| (see also El)- As was shown in [T^, the F{R) gravity can be treated as 
a perfect fluid, so the slow-roll parameters and the observational indices can be calculated in terms of the Hubble 
rate easily. In our case, the procedure is the same, so we do not reproduce it here. Hence, we shall calculate the 
slow-roll parameters, for a viable cosmological expansion and then determine the potential V{t = <f), which generates 
this expansion. In order to have a better command on the resulting equations, we shall use the e-folding number N^ 
instead of the cosmic time t. Using the following transformation rules for the derivatives with respect to the cosmic 
time and to the e-folding number. 




.i7(iV) —— 
^ ^ dA^ dA^ 


and following the slow-roll indices can be written in terms of the e-folding number N, as follows. 


(16) 


H{N) 

4:H'{N) 


( H"{N) 
H(N) 


V 


-b6 


H'jN) 

H(N) 



H'jN) 

HiN) 


2 


11 = - 


n H'(N) o H"(N) i ( H\N) Y _ i ( H" {N) Y q H" {N) H'" (N) 

H{N) H{N) 2 H{N) J 2 \ H'{N) J H'{N) H>{N) 


2(3 + 


H{N) ) 


(17) 


where the prime denotes differentiation with respect to the e-folding number N. The spectral index of primordial 
curvature perturbations and the scalar-to-tensor ratio r, are equal to. 


Ug ~ 1 — 6e-I-2ry, r = 16e . (18) 

Notice that these relations are valid, so long as during inflation (which means for most values of N), the slow-roll 
indexes are e,ri ^ 1, so this formalism is valid only if the slow-roll approximation is assumed to hold true. 
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Having these at hand, we can choose a cosmological evolution compatible with the recent Planck data and find 
the corresponding potential V{N) that generates this evolution, for the F{R) function being of the form (fT^ . As 
we already mentioned, we shall be interested in the spectral index of primordial curvature perturbations Us and the 
scalar-to-tensor ratio r, which are constrained by the recent Planck data [l^, as follows, 

= 0.9644 ±0.0049, r<0.10. (19) 


Let us now exemplify our method by using some characteristic examples. Consider for example the class of models of 
cosmological evolution with Hubble rate of the form, 

H{N) = {-Go + Gi)'’ , (20) 

where Gg, Gi, /3 and b arbitrary real numbers in general. Since for negative values of Gg and Gi, the Hubble rate 
might become negative, we choose b to be of the form, b = 2 ^^, with m and n positive integers chosen in such a 
way so that 6 < 1. The slow roll parameters for the Hubble rate (pH)) are equal to. 


e = N^-^(Gi-GoN^'^(s- 


bGoN-^+Pp 
Gi - GoNf^ 


) "(46Gg^)-i( 


QbGoN-^+Pp 
Gi - GoNP 


±(Gi-Go 7V^) ''(-&GgiV-2+/5(Gi-GoiV'5) '^'(-l±^)/3 
± (-1 ± h)hGlN-‘^+‘^^ (Gi - GoN^'^ ^ • 


62GgAf-2+2/3/32 

{Gi-GoNP)^ 


with regards to the parameter e, while the parameter rj reads. 


G?(-l ± /3)(-3 + m + P) + GgAf2/3 (3 _ io5;3 + 352^2 6 ^(_i + 4 ^^)^ 

AN (Gi - GoNf^) {3GiN - GoNl^{3N + b/3)) 
-2GgGiiV/5(3iV(-2 ± /3 ± 46/3) ± (-1 ± /3)(-3 ± (-1 ± 56)/3)) 

4Af(Gi-GoAf/3)(3GiA^-GoA^/5(3/v + 6/3)) 


In view of Eqs. dan and (I22L the observational indices of Eq. (IT^ become. 


( 21 ) 


( 22 ) 


Us = -2 ^ (23) 

2n(Gi - GqNP'^ ( - 3GiiV ± GoN0{3N + 6/3)) 

(5GIn(^ - 3 ± 6A^(1 + N)+AI3- 6Np - - GIN^^{9N{-1 + 2N{1 + N)) 

+ ASbN^P + 2b‘^{-l ± 13Af)^2 453^3^ _ ^54iv3 + 26(-l ± + 3A^(-1 ± ^)(9 ± /3) 

± 6N^{9 - 6/3 ± 86/3)) ± GIGiN"^!^(^5AN^ + 26(1 ± 6)(-l ± 

+ 6N^i9 - 3/3 ± 166/3) ± A^( - 27 ± 2/3(6 ± 3/3 ± 136^/3)))) . 


By choosing the free parameters Gg, Gi, /3 and 6 as follows. 


Go = 0.00005, Gi = 3000, /3 = 3, 6 = ^ , (24) 

9 

the corresponding spectral index and scalar-to-tensor ratio for nearly e-folds {N = 60), are equal to, 

n, ~ 0.966157, r ~ 0.000258947, (25) 

which are compatible with the observational constraints m- Hence, since the cosmology generated by (uni is 
compatible with observational data, we can straightforwardly find the scalar potential that can generate such a 
cosmology. By introducing the function G{N), which is equal to G{N) = H{N)'^, we can easily see that the Ricci 
scalar R, is related to the function G{N) as follows (see also 0), 


RiN) = 3G'{N) + 12GiN ). 


( 26 ) 
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By combining relations (l20l) . (|T5|) and (l26l) . and by substituting in (I13p . we obtain the scalar potential V{N), which 
is equal to, 

V{N) = 6 {IC{N)f'' + 12d + 432d {IC{N)f'' + 1728d (/C(iV))®'' (27) 

- 4bGoN-^+^ (/C(iV))"^+^''/3 - 2166dGoA^“^+'^ /3 - 2592bdGoN-^+^ (/C(iV))“^+®''/3 

+ 1296b^dGlN-^+^l^ (/C(A^))"^+“ (3^ - 216b^dGlN-^+^P (/C(7V))“^+®^/33 , 


with JC{N) = Gi — GqN^ and the corresponding Lagrange multiplier of Eq. (IT4)) . is equal to, 


\{N) = -3 (/C(iV))^'' - 216d {JC{N)f’’ + 432d (/C(iV))®^ 

+ 1086dGofV-i+^ (/C(iV))”^+'‘^^- 12966dGofV-i+^ (/C(Af))"^+“ ^ 

+ 972b^dGlN-^+^P (/C(A^))"^+“ f3^ - 216b^dGlN-^+^l^ (/C(A^))"^+“ l3^ . 


(28) 


Then, by recalling that (j) = t, we can in principle find the relation t = t{a), where a is the scale factor, and express 
the cosmic time as a function of N. However, in the present case this is not possible, since the function t = t{a) is of 
the following form. 


t = —F2 
^0 


/ 


1 

-''I ~n 5 


1 Go In 


{i)‘ 


V 




Gi 


In 


(29) 


which cannot be solved as a function of a, where F 2 is the Gauss hypergeometric function, but the potential ((27ll 
provides enough information for the behavior of the mimetic matter. There are cases however, that it is possible to 
find the potential as a function of (j>, for example for the cosmological evolution with Hubble rate, 


H{N)={-Goe^^ + Gi) , 

it is possible to find the potential as a function of (f). The slow-roll parameter e for this evolution reads, 

be^^Go^ (Gi(6 + ^)- 2e^^Go{3 + 6/3))" 


e = —- 


AFiN) 

with the function F{N), being equal to, 

F{N) = {e^^Go - Gi) (-3Gi + e^^Go(3 + 5^))" 
while the parameter 77 can be easily calculated and it is equal to, 

/3 (862e2/5^G2/3 + Gi (2e'5^Go(-3 + ;9) + Gi(6 + /3)) + 2be^^Go {Ue^^Go - Gi(12 + 5/3))) 


f] = 


4 {eP^Go - Gi) (-3Gi + e^^Go{3 + bP)) 


(30) 


(31) 


(32) 


(33) 


By combining Eqs. and (ESI), we can easily compute the corresponding spectral index and scalar-to-tensor ratio, 
with rig being equal to, 

2(e^)^'^G3(3 + 5^)2(l + 26/3)+3G?(-6 + 6/3 + ; 92 ) e^^^GoGf (54 + 12(-3 + 46)/3 + 3/3^ + 26/33) 

"• =- mm -+- mm - 

2e^/^^GlGi (27 + (-9 + 486)/3 + (3 + 136^) + 6(1 + b)P^) 

2F{N) ■ 

and the scalar-to-tensor ratio is equal to, 

Abe^^GoP (Gi(6 + /3) - 2e^^Go(3 + 6/3))" 


r = -- 


F{N) 


(35) 


By choosing the parameters Gq, Gi, /3 and 6 as follows. 


Go = 0.5, Gi = 10, P = 0.0222, 6 =1 


(36) 
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the indices Ug and r read, 


ug ~ 0.96567, r = 0.0640848, (37) 

which are in concordance with the Planck data m- Of course there are other values of the parameters Gq, Gi, /3 
and 6, for which the observational indices are compatible with observations, like for example if we choose. Go = 0.5, 
Gi = 12, /3 = 0.024 and b = 1/2, but for the choice (I36L it is possible to express the potential V{N), as a function of 
(j). Indeed, we can easily obtain the potential V{N), which is, 

V{N) = 6 + 12d + 3 + 27d (5(iV))^^'' + 26e'^^Go /3 (38) 

+ 432d + 1836d (5(Af))“ + 1296d {S{N))^'' + 324d (5(7V))^°'’ , 

with S{N) = —e^^Go + Gi, and the corresponding Lagrange multiplier reads, 

X{N) = -s(siN)y^ - 6d(s{N)y^ - - 18d(5(Af))^^ - 186de^'^Go(5(7V)) l3 (39) 

- ^d(5(iV)) - Sbe^^Go (s{N)^ - 72bde^^Go (s{N)^ ■ 

By choosing the parameters as in Eq. (1361) . it is possible to express the cosmic time as a function of the scale factor 
a, in the following way: since H{N{a)) is equal to, 

H{N{a)) = (-Go {a/aof + Gi)\ (40) 

owing to the fact that e'^ = ^, by integrating the above equation and for the values of the parameters given in Eq. 
(IMl) (5 = 1 is the most relevant for the calculation), we obtain the function t = t{a), which is, 


t + Cl — 


/31n(a)-ln(Gi-Go 


Gi/3 


which can be solved with respect to the scale factor a, 


y{t) = 


gGi/3(t+ci)^^ 

1 _|_ eGi/3(t+ci) Go 
“o 


Cl = 


■In 


Q!0 


/3Gi Gi - 


Recalling that ^ = t, by making the following substitution in Eq. 


odN 


gGi/3(0+ci)(^^ 

1 _|_ gGi/3(0+ci) ’ 


(41) 


(42) 


(43) 


we easily obtain the scalar potential In the context of the reconstruction method we used, it is possible to 

reproduce a large number of phenomenologically appealing cosmologies, by just using any viable F{R) gravity one 
may choose and also any given Hubble rate, which produces a cosmological evolution compatible with observational 
data. We used a simple F{R) function, but in principle one can choose an exponential F{R) gravity or other viable 
models 0. In this way, it is possible that the gravitational theory at hand, may have all the appealing attributes 
of F{R) theories, like explaining early and late-time acceleration, and at the same time providing an alternative 
candidate for dark matter. 


III. INVERSE RECONSTRUCTION METHOD FOR DETERMINING THE F{R) GRAVITY 

The reconstruction method we presented in this letter can work in the opposite way too, for example by providing 
the Hubble rate, the potential E(/)) and the Lagrange multiplier X{(j)), it is possible to find the F{R) gravity which 
generates such a cosmolo gica l evolution. Let us exemplify this by using a symmetric bounce, which is a very well 
known bounce cosmology [l^, with the scale factor and the corresponding Hubble rate being equal to, 

a{t) = e'*'* , H(t) = 2"/t. 


(44) 
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Suppose that the potential V{4> = t) and the Lagrange multiplier X{t) are given as follows, 

F(t) = -24tV, A(t) = -27. 


(45) 


Having at hand the cosmological evolution (l44ll . the potential and the Lagrange multiplier of Eq. (145 1) . by making 
use of Eqs. (HU and da, we may easily obtain the F[R) gravity that generates such a cosmological scenario. Upon 
combining Eqs. 0 and da, we obtain the following differential equation. 


d^F'{R) 


dF'{R) 


AF'{R) 


dF'{R) 


2 " ~ + 477 -^^)"^ - 6H—^ - 2{H + + 6(H + _ 2X{t) - V{t) = 0 (46) 

at^ at at at 

Using Eqs. and ([35]), and applying these to Eq. ([i5|) . we obtain the following differential equation, 

2y{t) + AH{t)y{t) - QH{t)y{t) + 6 (H{t) + Hitf) y{t) - 2 - 2X{t) - V{t) = 0 , (47) 

where y{t) = F'{R{t)), and solving this differential equation with respect to y{t), one gets the following solution. 


yit) = Cie-K5+V37)*^^^^ (371/4^^) ^ ^ i, y^ty), 


(48) 


where ni = ^ (—37— Cl and C 2 arbitrary constants, and the functions Hn{x) and iEi(a, b, x) are the Hermite 

polynomial and the Kummer confluent hypergeometric function respectively. By also recalling that R = 677 + 1277^, 
substituting the Hubble rate (l44l) and solving with respect to t, we obtain, t = , so by substituting this in 

(l48l) . we can easily obtain the function y{R). Then we can obtain the analytic form of F{R) in the large and small R 
limits. For example, for large 7?, the F'(R) gravity is approximately equal to, 

7^'(7?)~Ci+C27?, (49) 


where Ci and C 2 are constants, the analytic form of which is irrelevant to our analysis and can be found in the 
Appendix A. 


IV. EXISTENCE OF DE SITTER SOLUTIONS AND DYNAMICAL EVOLUTION 

The early time acceleration era is described by (quasi)-de Sitter solutions, so firstly it is important to investigate if 
de Sitter solutions exist in the case of mimetic F{R) theory with potential V{t = (f)). If a specific F{R) has de Sitter 
solutions, we need to examine whether these de Sitter solutions are stable against linear perturbations. If the de Sitter 
solutions are stable, then this solution is the final attractor of the theory and therefore the eternal inflationary era 
lasts forever [T^ . and on the contrary, if the perturbation is unstable, the graceful exit from inflation is achieved. 
We start off with the F{R) gravity of Eq. (fT51) . For this model, de Sitter solutions exist only for d < 0, as we now 
demonstrate. The potential V{t) and the Lagrange multiplier X{t), for the Hubble rate being 77(t) = H^s, take the 
following form, 

V{t) = mis - X{t) = -mis - 432d776g , (50) 

and by substituting these to Eq. (1461) . we get the equation QFlls (—1 + 432d77jg) = 0, and by solving it, with respect 
to 77ds, we obtain the only de Sitter solution that exists, which is of the form, 

= 2 33/4(71/4 ■ 

What now remains is to study the linear stability of this solution. An unstable de Sitter solution indicates that the 
inflation ends. We consider perturbative solutions of the de Sitter solution H{t) = Flj^St which have the following 
form, 

H{t) = Hds + A77(t), (52) 

with the perturbation being very small, that is | A77(t) |«C 1. Substituting this to Eq. (HSl) . and by keeping terms 
linear to AH{t) and it’s higher derivatives AH{t) and A77(t), we obtain the following differential equation for the 
perturbation AH{t), 

-2592dHlsAH{t) - 12HdsAH{t) + 15552d77j5A77(t) - 6A77(t) + 5184d77^5A77(t) (53) 

-10368d77JsA77(f) - 677Js + 2592d77^s = 0 









where the “dot” here denotes differentiation with respect to the cosmic time t. Note that we used again Eqs. (nsi) 
and (HI to which we substituted Eq. dSll). The solution to this equation is of the form, 

AH(t)=C+Cie‘^^*+C2e^^\ (54) 

where Ci,C2 are arbitrary integration constants and the exact analytic form of the parameters C, Ci and C2 can be 
found in the Appendix, since these are too complicated to be presented here. By substituting the exact de Sitter 
value Hds of Eq. dSB), in the parameter C 2 , we can see that the parameter (^2 is positive, so the term ~ dominates 
the evolution, since the other term contains (( 1 , which is negative. Therefore, the de Sitter solution we found in the 
F{R) = R — dR? model is unstable, and therefore inflation can end in this model, always in the presence of mimetic 
potential and Lagrange multiplier. 

Before closing our study, it is worth presenting the results of another model, for which a de Sitter solution exists 
and the perturbation evolution is much more unstable in comparison to the model we just studied. The model is 
described by the following F{R) function, 


F{R) =R-dR^ + fR^ , (55) 

so it results from the model of Eq. dig, by simply adding an R^ term. Interestingly enough, the model (1551) is known 
to yield quite interesting solutions for compact objects as was shown in [^ . The intriguing feature is that when the 
F{R) gravity is taken exactly as in Eq. (1551) (that is, the coefficient of R^ negative and the coefficient of R^ positive), 
there appears an increase in the maximal neutron star mass. The full study of the model (1551) will be presented 
elsewhere, here we describe only the results of our analysis. The de Sitter solution for the model (|55|) is the following, 

„ !,/./ , VsdTp 

d ’ 

and the evolution of linear perturbations behaves as, 

AH{t) =A + CAe^^^+ CBe^^\ (57) 

where A, 61 and S 2 appear in the Appendix, and the parameters ca and Cb are integration constants. The parameter 
i5i is negative, and the parameter S 2 is positive. In addition, S 2 grows larger as / and d appearing in Eq. (1551) increase. 
The resulting picture is that the linear perturbations of the de Sitter solution of the model (l5^ are much more 
unstable in comparison to the model dig, which means that the R^ term makes the model more unstable towards 
linear perturbations of the de Sitter solutions as it was demonstrated long ago by Starobinsky [2l|. The detailed 
study of this model will be given elsewhere. Finally, we need to mention that the model F(R) = R + dR? + fR^ has 
no de Sitter solutions, so only the model with negative d and it’s R^ deformations have de Sitter solutions. 


V. CONCLUSION 

In this letter we studied the behavior of mimetic F{R) gravity with scalar potential. As we demonstrated it 
is possible to generate inflationary cosmology consistent with observations by using specific forms of the mimetic 
potential and Lagrange multiplier, for a given form of the Hubble rate. After exemplifying our findings by using some 
characteristic examples, we investigated the existence of de Sitter solutions for the mimetic potential F{R) theory and 
also we examined the stability of these de Sitter solutions. As we evinced, the final result strongly depends on the 
form of the F{R) gravity. It is worth to study the existence of de Sitter solutions for specific classes of potentials. This 
procedure would narrow down the possible F{R) gravities that would yield de Sitter solutions for the given potential. 

In addition, the Einstein frame counterpart theory of the mimetic F{R) theory can easily be found, and the final 
form of the resulting two scalar theory is interesting to study. Particularly, since in the Einstein frame there is a direct 
coupling between the Einstein frame scalar and the mimetic scalar, it would be interesting to investigate the resulting 
theory, which can have similarities with some recently studied two scalar field models [ 2 ^ , if the Einstein frame scalar 
counterpart of the F{R) gravity is chosen to be the R^ inflation model Also the issue of gravitational particle 

production should be re-addressed in such a theoretical framework. 

Moreover, it would be very interesting to reproduce known consistent with observations Jordan frame cosmologies, 
such as the R^ inflation cosmology, by using an arbitrary viable F{R) gravity and a suitably chosen mimetic potential. 
Finally, a study of the dynamical system that corresponds to the F{R) gravity with mimetic potential would be of 
importance. Specifically, it is worth investigating which are the fixed points of the theory, also perform a stability 
analysis of these points, and also providing the physical meaning of some of these fixed points could be of importance. 
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since the mimetic potential could have interesting implications on the physical meaning of the fixed points. Finally, 
we need to note that in this paper we assumed a kind of slow-roll approximation in order to extract the observational 
indices in the Jordan frame. It is p ossible however to abandon the slow-roll approximation and to study different 
inflation regimes (compare with |24|1. 

Finally, we need to note that the mimetic F{R) formalism can be applied in such a way so contact with the present 
Universe is achieved. Particularly, this can be done by using some recent cosmographic developments [Ull^. There 
are three possible ways to accommodate cosmographic data in the mimetic F{R) formalism. Firstly, one could use 
the resulting Hubble rate of cosmographic data, which is of the form [2^, 

H{z) ^ + ln(a -I- /3z) (58) 

where 2 ; is the redshift, and a and /3 parameters of the theory (see [1^ for details), and for a general F{R) gravity it 
is possible to find the Lagrange multiplier X{t) and the potential V{t) that can yield the evolution (1551) . The second 
way that contact with cosmography can be achieved, is to use the F{R) gravity which corresponds to the latest 
cosmographic data [ 2 ^ (see Eq. 31 of [1^), and find the potential and Lagrange multiplier that generates such a 
cosmological evolution. Finally, the third way to make contact with the latest cosmographic data is to use both the 
Hubble rate and the F{R) gravity that correspond to cosmographic data, and again investigate how the potential 
and the Lagrange multiplier behave for these choices. Lastly, let us note that it is quite interesting to generate the 
cosmological evolution of Eq. ((55)) by using the mimetic F{R) gravity formalism, with the F{R) gravity being slightly 
different from the General Relativity case, for example F{R) ^ with e <C 1. In this way we could find a power 

expansion of the Lagrange multiplier and of the potential, in terms of the parameter e and the resulting expressions 
would be quite simple. We hope to address these issues in a future work. 
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Appendix: Detailed Expressions of Parameters Appearing in Text 


Here we quote the exact form of the parameters Ci, C 2 , C, ^i, ^ 2 ) -4, di and <52 appearing in the main text of the 
paper. We start off with Ci and C 2 which are the coefficients of the small curvature expansion of the F'{R) gravity, 
appearing in Eq. (l49l) . the analytic form of which is, 


Cl = 




(59) 


C 2 = Cl 





1 277 ) 

1 (37+ Iiy^) I 

(^) 

l\ 

1 967 


48 373 / 473/2 



'(5 + v^)e"i 


C 2 


7i^^i (37 + e-3(5+'^)^iFi (1 - 


967 


96v^7 


where Ua is equal to, 


na = ^(-37-IlV37) . 


(60) 
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The parameters C, Ci and C 2 that appear in Eq. (I54L are equal to, 


C = 


Ci = 


C2 = 


Has (-1 + 432di?|5) 

2 (1- l2mdHjg) ’ 

(-1 + 864di7jg - n2MHlg - - 1728dH^s + 746496d2Ef«s + 1492992d2Efys + 2985984d2Efi^) 

86Wj^ 

(-1 + S64dHjs - l728dH^s + y'l - l728dH^s + 746496d2Er«s + 1492992d2Er«s + 2985984d2Efj^) 

86 Wj^ 


(61) 


Finally, the parameters A, Si and S 2 appearing in Eq. dSZJ, are equal to, 


, ^ (-l-24dgjg + 432dg4^) 

2(l + 48difJ5- 1296dij4^) ’ ^ ^ 

<5i = „ , ^ ,rr4 X f - 1 - 36dHL + 48dHL + 864dHL - 1728dHL 

2 (-12dH^g + 432dff;^g) \ 

- \J-4 (-12diJ2^ + 432dHjg) {2Hds + - 2592dH^g) + (1 + 36dHjg - 48dHjg - 864dHjg + 1728dH^g f'^ , 

S 2 = „ , ^ , f - 1 - 86dHjg + 48dHjg + 864dHjg - 1728diJ5c 

+ a /-4 {-12dHjg + 432dHjg) {2Hds + 9QdH^g - 2592dH^g) + (1 + 36dH^g - 48dH^g - 864dHjg + 1728dH^gf^ . 
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